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Mean-field limits for some Riesz interaction gradient flows 


Mitia Duerinckx 


This paper is concerned with the mean-field limit for the gradient flow evolution of particle systems with 
pairwise Riesz interactions, as the number of particles tends to infinity. Based on a modulated energy method, 
using regularity and stability properties of the limiting equation, as inspired by the work of Serfaty [28] in the 
context of the Ginzburg-Landau vortices, we prove a mean-field limit result in dimensions 1 and 2 in cases for 
which this problem was still open. 


1 Introduction 


We consider the energy of a system of N particles in the Euclidean space (d > 1) interacting via (repulsive) 
Riesz pairwise interactions: 


N 


HN{xi,...,XN) = '^9s{Xi-Xj), Xi€. 


where the interaction kernel is given by 


9six) := 


c^ jjxj if 0 < s < d; 

,-Cd,dog(|a^l)^ ifs = 0; 


with Cd,s > 0 some normalization constants. We note that the Coulomb case corresponds to the choice s = d — 2, 
d > 2. Particle systems with more general Riesz interactions as considered here are extensively motivated in 
the physics literature (cf. for instance [2, 20]), as well as in the context of approximation theory with the study 
of Fekete points (cf. [14] and the references therein). Recently, a detailed description of such systems beyond 
the mean-field limit in the static case was obtained in [21], and also in [17] for non-zero temperature. In the 
present contribution, we are rather interested in the dynamics of such systems, and more precisely in a rigorous 
justification of the mean-field limit of their gradient flow evolution as the number N of particles tends to infinity, 
which has indeed remained an open problem whenever s>d—2,s>0,d>2. 

We thus consider the trajectories x* ^ driven by the corresponding flow, i.e. the solutions to the following 
system of ODEs: 

9txlj^ =——'^iHiy{x\j^,... ,x%pf)^ xl^j^\t=o = x°^]^, i = (1-1) 


where (a;° 7 v)fci ^ sequence of N distinct initial positions. Since energy can only decrease in time and since 
the interaction is repulsive, particles cannot collide, and moreover it is easily seen that a particle cannot escape 
to infinity in finite time; from these observations and from the Picard-Lindelof theorem, we may conclude that 
the trajectories xl are smooth and well-defined on the whole of M+ := [0, cx)). As the number of particles gets 
large, we would naturally like to pass to a continuum description of the system, in terms of the particle density 
distribution. For that purpose, we define the empirical measure associated with the point-vortex dynamics: 



i=l 


(1.2) 


1 


and the question is then to understand the limit of as f oo. More precisely, assuming convergence at initial 
time 

1 ^ 

]V ^ as IV too, 

i=l 

formal computations lead us to expect under fairly general assumptions /r* for all t > 0, where /i* is a 

solution to the following nonlocal nonlinear PDE on M’*' x sometimes called the fractional porous medium 
equation: 


dtfJ* = div(/r‘Vh*), h* := gs * M‘|t=o = M°- (1-3) 

This equation in the weak sense just means the following: /i € lL(k+;l 1(M‘')), g. */i e lL(k+;<oc (»"))> 
/iVgs * g € Lj'Q(,(R+; L^(]R‘^)), and, for all </> £ (^“(IR^)) such that cj)(t, •) = 0 for alH > 0 large enough, 

/ / fj*{x){dt(l){t,x) —'V(f>{t,x) ■ ygs * f/{x))dxdt + / fj°{x)(l){0,x)dx = 0. 

Js.+ JR'i JR'i 

As far as existence issues as well as basic properties of the solutions of (1.3) are concerned, we refer to [9, 8] for 
d — 2 < s < d, s > 0, to [19, 12, 1, 29] for s = d — 2, d > 2, and to [10] for 0 < s < d — 2. See also Proposition 2.3 
below. 

In the case of logarithmic interactions s = 0, this expected mean-field limit result was essentially first proven 
(in arbitrary dimension) by Schochet [27[ based on his simplification [26] of the proof of Delort’s theorem [11] 
on existence of weak solutions to the 2D Euler equation with initial nonnegative vorticity in Schochet’s 

original paper [27] was actually only concerned with the mean-field limit for a particle approximation of the 2D 
Euler equation, but the same argument directly applies to the present setting. However, due to a possible lack 
of uniqueness of weak solutions to equation (1.3), Schochet [27] could only prove that the empirical measure 
/i(y converges up to a subsequence to some solution of (1.3). The key idea, which only holds for logarithmic 
interactions, consists in exploiting some logarithmic gain of integrability to find uniform bounds on the number of 
close particles, which allows to directly pass to the limit in the equation and conclude by a compactness argument. 

In the case 0<s<d — 2,d>3, the complete mean-field limit result (not restricted to a subsequence) was 
proven more recently by Hauray [15] (see also [10]), but his method, based on a control of the infinite Wasserstein 
distance, cannot be adapted at all to higher powers s > d — 2. In the ID case, Berman and Onnheim [3] obtained 
a similar result for the whole range 0 < s < 1, in the framework of Wasserstein gradient flows, but their method 
cannot be extended to higher dimensions. 

Very recently, in the context of the 2D Gross-Pitaevskii and parabolic Ginzburg-Landau equations, Serfaty [28] 
proposed a new way of proving such mean-field limits ^ , based on a Gronwall argument for the so-called modulated 
energy, which is some adapted measure of the distance to the (postulated) limit. This idea of proof originates in 
the relative entropy method first introduced by Yau [31] for hydrodynamic limits (see e.g. [23] and the references 
therein for later developments), the modulated Hamiltonian method used by Grenier [13] for boundary layer 
problems, and the modulated energy method designed by Brenier [6] for the quasi-neutral limit of the Vlasov- 
Poisson system. The advantage of this method is to be completely global, bypassing the need for a precise 
understanding of the microscopic dynamics. It relies on the regularity of the solution to the limiting equation, 
and exploits its stability properties. As will be seen, however, we are able to apply this method in the present 
context only in dimensions 1 and 2 and for s not too large. More precisely, we treat in ID the whole range 
0 < s < 1 as in [3], while in 2D we treat but the case 0 < s < 1, which is new and in particular completes 
Schochet’s partial result [27] in the logarithmic case. Our main result, for which we need an additional regularity 
assumption on the limiting equation (cf. Remark 1.2(a) below), is as follows: 

Theorem 1.1. Let d = 1 or 2, and 0 < s < 1. Let be a probability measure, and in the case s = 0 also assume 
Jgrf log(2-|-|a;|)|/r°(a;)[(ia; < oo. Assume that equation (1.3) admits a solution jjf that belongs to L°°([0, T]; C"^(]R'^)) 
for some T > 0 and some a > 2 — d + s. In the case s = 0, d = 1, also assume V/r £ L°°([0, T]; L^(R'^)) for some 

1. In [28|, the questions are different in nature, since they consist in passing to the limit in PDE evolutions, but are similar in spirit 
since one wishes to understand the limiting dynamics of point vortices which essentially behave like Coulomb-interacting particles. 
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(1.4) 


p < oo. Let p,°j^ ^ p° as above, assume the eonvergence of the initial energy 

.. ,x% m) = I [ gs{x - y)dfi°{x)dp°{y) < oo, 

^too ’ 7R<i jRrf 

and let be defined by (1.1)-(1.2). Then pf is the only weak solution to (1.3) up to timeT m L°“([0, T]; L°°(R'^)), 
and for all t € [0,T] we have ^ pf, as well as the convergence of the energy 

■ ■ ■>^n,n) = I 9s{x - y)dp\x)dp\y) < oo. 

^too iV^ JKd- Jm<1 


Remarks 1.2. 

(a) For a compactly supported probability measure p° S L'^(]R^), equation (1.3) always admits a solution in 
L°°(R+; L°°(R‘^)), which remains a compactly supported probability measure for all times (see Proposition 2.3 
below). As far as the additional regularity assumption is concerned, as explained in Section 2.3, it has been 
proven to hold with T = oo in the case 0<s<d — 2, d>2, and at least up to some time T > 0 in the case 
d — 2 < s < d — 1, s > 0, for sufficiently smooth initial data p°, but all other cases remain unsolved, and the 
additional regularity assumption in Theorem 1.1 above is then crucially needed. 

In the 2D Coulomb case s = 0, d = 2, as this regularity problem is solved (cf. [19, Theorem 1]), the conclusion 
of Theorem 1.1 above holds in that case with T = oo (that is, with L°“([0,T]; ■) replaced by L“(,(]R+; •)), 
whenever we have p° £ 7^(]R^) n C'^(M^) for some a > 0, and log(2 + |a;|)|/i°(a;)|da; < oo. This completes 
Schochet’s partial result [27]. As the regularity problem is further solved in short time in Sobolev spaces in 
the case d—2<s<d — l,s>0 (cf. [30[), the conclusion of Theorem 1.1 above holds for some T > 0 
(depending on initial data) in the case 0<s<l,d = 2, whenever p° £ P(]R^) fl for some cr > 2, 

and also in the ID logarithmic case s = 0, d = 1, whenever p° £ 7^(R) n id°'(R) for some a > 3/2, and 
/Rlog(2 + |a;|)[/r°(x)|da: < oo. 

(b) A closer look at the proof actually shows the following quantitative statement, where the distance between 

and p* is measured in terms of the modulated energy: for all t £ [0, T], we have for some Ct > 0 


< 


// 9six - y)d{p% - p*)ix)d{p% - p*){y) 

JJ x^y 

Ct [[ 9s{x- y)d{p% - p°){x)d{p% - p°){y) + Ct 

JJ x^y 


N l + s-o- 
N-^\og N, 


if s > 0; 
if s = 0. 


(c) The well-preparedness assumption (1.4) for the initial positions {x° N >0, is statistically relevant, in 

the sense that it is automatically satisfied almost surely if, for each N, the positions {x° are chosen 

to be independent R'^-valued random variables with law p°. This easily follows from the strong law of large 
numbers, together with the bound (for s > 0) 



y)\dp°{x)dp°{y) < 


|x-y|<l 


\x-y\ dp°{x)dp° {y) 


\x-y\>l 


dp°{x)dp°{y) < ||Ai°[|L“ 


+ 1 . 


(d) We may also add to the energy a potential V, thus turning Hjsr into := Hjsr + V. li V £ C'^(M'^) satisfies 
||V^D||l“ < oo, then all the arguments may be directly adapted, as long as the corresponding limit equation 


dtp* = div{p*V{h*+ V)), h*:=gs*p*, p*\t=o=p°, 


admits a regular enough solution. 

(e) In dimension d = 2, we could also consider a mix between the gradient flow (1.1) and its conservative 
counterpart, that is, replacing (1.1) by the following system of DDEs, for i = 1,..., A^, 


dtx\ jq — ——'^iHN{x\^i\[, ■ 


Xn,n) ~ 


t 

l.ATi ■ 


jXn.n) ~ m). 


: I _ o 

;jv|t =0 — V; 


3 



where we have also added a potential V as in item (d) above. If a > 0, then all the arguments may again be 
directly adapted, as long as the corresponding limit equation 

dtpi^ = div(fj,^(aVh* +/3V^h*+VV)), h* := /r‘|t=o = 

admits a regular enough solution. (Note that the same proof can a priori not work for the choice a = 0, since 
in Step 2 of the proof of Proposition 2.8 below some term cannot be estimated directly and needs instead to 

be absorbed using the negative diffusion term, which would however vanish in the case a = 0.) 

Notation. Denote by B{x,r) the ball of radius r centered at x in and set Br := B{0,r). We also use the 

notation a A b = min{a, b} and aV b = max{a, b} for all a,b € M.. The space of probability measures on is 

denoted by 'P(R'*), and, for all tr > 0, C'°'(R'^) stands as usual for the Holder space (R'^), while C'(f(R'^) 

denotes the subspace of compactly supported functions. In the sequel, C denotes any positive constant only 
depending on d and s. We denote by Ct any positive constant only depending on d, s and on time t, such that 
Ct < Ct for all t G [0, T] and all T > 0. We also denote x < y and x <t y ioi x < Cy and x < Cty, respectively, 
and we use the notation x ~ y if both x < y and y ^ x hold. Finally, we denote by Oa(l) a quantity that goes 
to 0 when the parameter a goes to its limit, uniformly with respect to other parameters, and we write oi^^(l) if 
it converges to 0 only for any fixed value of the parameter b. 


2 Proof of Theorem 1.1 

2.1 Strategy of the proof 

Translating the idea of [28] in the present setting (see also [6]), the clue of the proof of Theorem 1.1 comes 
from the following stability estimate, which we first present for simplicity in the Coulomb case. 

Lemma 2.1 (Stability — Coulomb case). Let s = d — 2, d > 2. Let £ T’{W’‘) 0 L°°(R'^), and in the 

case d = 2 also assume log(2 + |x|)(|p,)’(x)| + \y 2 {x)\)dx < oo. For z = 1,2, let yl be a weak solution of 

equation (1.3) with initial condition y,°, denote h\ := gd -2 * lAi o-nd assume /ri,/r 2 £ L°“([0, Tj; L°“(R‘^)) and 

V^/i 2 £ L^([0, T]; L°°(M‘^)) for some T > 0. Then, for all t £ [0,T], 

[ [ gd- 2 [x-y)d{y\-y\){x)d{y\-yl){y) (2.1) 

< WV^KW^^du f f 
iR'i JR‘‘ 

Proof. Proposition 2.3(ii) below yields V(/ii — / 12 ) £ L°°(R+; L^(R'^)). Combining this with the additional bound¬ 
edness assumptions, all integration by parts arguments in the sequel may be justified. 

Using the equations for and the time derivative of the left-hand side of (2.1) can be computed as follows: 

dt / 5d-2(a: - y)d{yl - y\){x)d{y\ - y\){y) 

= 2 [ {h[ - hl){dty{ - dtyl) = -2 [ V{h[ - hl){fi\Vh{ - y* V/z‘) 
iR'i Jr<‘ 

= - 2 [ \Vihl - /z‘)|V‘i -2 [ Vhl- V(/z‘ - hi) (y‘ - yl). (2.2) 

JR'l dR<i 

The first term in the right-hand side is nonpositive, so it suffices to estimate the second one. Using the relations 
-Ah* = yl, i = 1,2 (which hold with a unit factor for the suitable choice of the normalization constant Cd,d -2 > 0), 
the product V(h* — hi) {y\ — yl) may be rewritten a la Delort using the stress-energy tensor: 

-2V(h*i - hi) {y\ - yl) = 2V(h* - h*) A(h‘ - h*) 

= div (2V(h* - h‘) 0 V(h* - hi) - Id |V(h* - h‘)p) , (2.3) 
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where the divergence of a 2-tensor here denotes the vector whose coordinates are the divergences of the corre¬ 
sponding columns of the tensor. Combining this with an integration by parts, we find 


2 [ • V(/i‘ - hi) (/r‘ -fj.l) = - f div (2V(/i‘ - hi) ® V(/i‘ - hi) - Id |V(/i‘ - • Vh: 

= [ (2V(/i‘i - hi) (g) - hi) - Id - hl)\^) : 

The inequality 2|a6| < -f and an integration by parts then yield 

f vhl ■ y{h\ - hi) - fil) < Wv^hlWi^^ [ \v{h\ - hl)\^ 

Jr<‘ Jr<‘ 

= I|W2IIl<» [ {h\-hl){fi\-f,l) 

Jr'^ 

= ||V 2 h* 2 ||L“ / / gd- 2 {x - y)d{g\ - yl){x)d(y\ - yl)[y), 

Jr‘^ 

so that the result (2.1) follows from (2.2) and a Gronwall argument. 


(2.4) 

□ 


We are also interested here in the non-Coulomb case d—2<s<d, s>0, and hence, just as in [21], we 
need to use the extension method of Caffarelli and Silvestre [7] (cf. Section 2.2 below) in order to find a similar 
Delort-type formula as in (2.3) of the proof above, and then repeat the same integration by parts argument, thus 
circumventing the fact that the Riesz kernel is not the convolution kernel of a local operator. This allows to prove 
the same estimate as above in all cases 0V(d — 2)<s<d with gd -2 replaced by g^ (cf. Lemma 2.2 below). 

This stability estimate gives us a control of the i7“^-distance (or the iJ“(‘^“'*)/^-distance, for general 0 < s < d) 
between fi\ and in terms of the initial distance, up to a factor that only depends on the regularity of in 
the form of ||V^h|||L°°. We would then like to replace by the smooth solution p* and to replace p* by /i)^. 

However, the corresponding distance would then be infinite because of the presence of Dirac masses, and moreover 
does not exactly satisfy the limiting equation (1.3). The idea of the proof of Theorem 1.1 consists in finding 
a suitable way of adapting the proof above to that setting. 

First, the natural way of giving a meaning to this divergent distance between and yd simply consists in 
excluding the diagonal terms, thus considering the (renormalized) modulated energy 


^Nit) = // gs{x - y)d{y% - y^){x)d{y% - p*)(y). 


I x^y 

The goal is then to compute the time-derivative dt^Nit), and trying to adapt the proof of the stability estimate 
above to bound it by CSnit) for some constant C > 0, up to a vanishing additive error. However, at the end 
of the proof above, the use of the inequality 2|a6| < is clearly not compatible with the removal of the 

diagonal terms. To solve this main issue, the crucial idea is due to Serfaty [28]: regularizing the Dirac masses 
at a (small) scale rj so that the diagonal terms become well-defined and diverge only as ry 4- 0, we need to try 
to construct around the particle locations small balls that contain most of the divergent ry-approximate energy, 
so that excluding diagonal terms essentially amounts to restricting the 7y-approximate integrals to outside these 
small balls. Using the same approximation argument as in [28] to be allowed to restrict all integrals to outside 
these balls, the end of the proof above is then easily adapted, using the inequality 2\ab\ < only on the 

restricted domain. 

In this way, for any 0V {d—2) < s < d, we manage to prove dt£N{t) < C£N{t) + on{^) under some mesoscopic 
regularity assumption on the distribution of the particles in time (cf. Proposition 2.8 below). Finally, in the case 
s < 1 (hence our limitation to that regime), these conditions can be directly checked using a modification of the 
ball construction introduced by [24, 16] for the analysis of the Ginzburg-Landau vortices (cf. Section 2.7 below). 

2.2 Extension representation for fractional Laplacian 

We recall here the extension representation for the fractional Laplacian by Caffarelli and Silvestre [7[ (we 
follow notation of [21, Section 1.2[). Let 0 V (d — 2) < s < d be fixed. For a finite Borel measure y on the 
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associated Riesz potential (—A) can be written (up to a constant) as 


h>^{x) := [ \x-z\ '‘dfj.iz) = gs * nix). (2.5) 

’ JR’i 

We denote coordinates in x R by (x,^), and we denote by n^R^i-xio} the Borel measure on R'^ x R defined as 
follows: for all (j) G C)(°(R‘^ x R), 


/ 4'{x,Odin5R<i^{Q})ix,C)-= <(>(x,0)d^(x). 

Jr-^xr Jr'^ 


Extending to R'^ x R via 

/i^(x,^):= f \ix,0 - iz,0)\~'^dn{z) = gs* in5Rd^{oy){x,0^ 

Jr‘>‘ 

where we denote gs{x,^) = c))'^|(x, 5)|“®, and choosing 7 :=s + l — ds (—1,1), the extended function on 
R^ X R satisfies in the distributional sense 


-div(|^PV/i'") = Ai4<Jx{o}- 

The function gs is indeed a fundamental solution of the operator — div(|^|'’'V) on R'^ x R, in the sense that 

— div(|^|^V(?s) = do on R'^ x R. The normalization constant Cd,s is chosen exactly to satisfy this property with a 
unit factor (for an explicit value, see Step 1 of the proof of Lemma 2.11 below). 

In the case s = 0, d = 1, denoting go{x,^) = —log(|(x,^)|), we have —Ago = on the extended space 
R X R, for the suitable choice of the normalization constant ci o, so the above again holds with 7 = 0 = s + l — d. 
(In the Coulomb case s = d — 2, d > 2, no extension is needed, and the normalization Cd,d -2 is simply chosen 
such that —Agd -2 = do on R'^.) 

Using this extension representation, we may now directly adapt the stability estimate of Lemma 2.1 to the 
non-Coulomb case: 

Lemma 2.2 (Stability — Riesz case). Let 0V(d —2) < s < d. Let/i 2 G 7^(R‘^)nL'^(R‘^), and in the case s = 0 
also assume log(2-|- |x|)(|/r)’(x)| -|- \n 2 ix)\)dx < oo. For t = 1, 2, let n\ be a weak solution of equation (1.3) with 
initial condition fj,°, denote h\:=gs*n\, and assume /xi,/X 2 G L°°([0, T]; L°°(R‘^)) and V^ti 2 G L^([0, T]; L°°(R‘^)) 
for some T > 0. Then, for all t G [0, T], 

[ f 9six - y)din\ - nl)ix)din\ - nl)iy) (2.6) 

Jr‘>- Jr’I' 

< [ [ g,(x - y)dinl - n°,)ix)dinl - g^)(g). 

dR<l jRd- 

Proof. By Lemma 2.1, we only need to consider the case d — 2<s<d, s>0. Proposition 2.3(ii) below yields 
V(tii — ti 2 ) G L°°(R+; L^(R'^ x R, |^|^dxd^)). Combining this with the boundedness assumptions, all integration 
by parts arguments in the sequel may be justified. 

Just as in (2.2), the time derivative of the left-hand side of (2.6) is 

dt [ [ gsix-y)d{n\-nl)ix)din\-nl)iy) (2.7) 

jR'i JR’‘ 

= - 2 [ \v{h\ - d‘)|V‘i - 2 [ vh\ • v(d‘ - h\) in{ - nl). 

Jr-^ Jr’I' 

The first term in the right-hand side is nonpositive, so it suffices to estimate the second one. Using the relations 

— div(|^|'''V/i‘) = /x‘dRdx{o}) foi' * = li2, we find the following proxy for the Delort-type formula (2.3): for all 
1 < fc < d, 

-2dk{h\ - hi) (/XidR<ix{o} “ d24<ix{o}) = 25fc(M - 4) div(|^pV(4 - hi)) 

d+1 

= (2|erafe(M - hDdlihi - hi) - 4/ieriv(4 - d‘)|2). 

1^1 
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Combining this with an integration by parts, we obtain 

« d d-\-l p 

2 / Vft.2 ■ V(/i* — / 12 ) (/r* — ^\) = — EEL di {2\i\^dk{h\ - h\)di{h\ - h\) - Ski\^r\v{h\ - hi)\^) dkh\ 

fc=l 

d d-\-l p 

L ler {2dk{hl - hl)di{h\ - hi) - 6ki\V{hl - hl)\^) dkihl 

k=i 1=1 •'R'^xr 

Hence, arguing as in Lemma 2.1, an integration by parts yields 


[ Vhl ■ V{hl - hi) (pI - pD 

< l|V2/t‘||L~ 

[ \^P\V{hl-hl)\^ 


dR<i 


JR’^xR 



= ||V2/t‘||L- 

/ (^1 ~ ^D(/^i'^R‘^x{ 0} “ /f2^R‘^x{0}) 




dR<ixR 



= wv^hlu^ 

f {h\ — hl){p\ — pI) 




dRd 



= l|V2/t‘||L~ 

/ / 9s{x-y)d{p\-pl){x)d{p\-pl){y), 

dRrf dRd 

(2.8) 


and the result (2.6) follows. □ 


2.3 Properties of the fractional porous medium equation 

Let us first state for weak solutions to (1.3) an existence result as well as some basic properties. We refer 
to [9, 8] for d — 2 < s < d, s > 0, and to [19, 1, 29, 4] for s = d — 2, d > 2. For 0 < s < d — 2, we refer to [10, 
Section 4], but there existence is proven only for short times, while in the present repulsive context it can easily 
be extended to all times, using that no blow-up can occur in finite time. ^ 

Proposition 2.3 (Existence for the fractional porous medium equation). Let 0 < s < d. 

(i) Existence: Let S P(]R‘^)nL°“(]R‘^), and in the case d—2 < s < d, s > 0 also assume that \fi°{x)\ < 

for some a,A>0. Then, there exists a (global) weak solution fT to (1.3) in L°“(]R"''; P(R‘^) n L°°(IR‘^)), 
which is unique in this class in the case 0<s<d — 2,d>2. 

(ii) General properties: Let p° G P(]R‘^)nL°°(]R‘^), and in the case s = 0 also assume log(2-|-|a;|)|/r°(x)| < c». 
Any weak solution pf to (1.3) on K’*' X with initial condition p° then satisfies 

/ / gs{x-y)dp\x)dp\y)< / gs{x - y)dp°{x)dp°(y), 

dR-i Js.'i Js.d Js.'i' 

where the left-hand side remains finite. Moreover, for all t > 0, 



y)d{p^ 


p°){x)d{p* 


T°){y) 


Jg,\V{h*-h°)\^, ifs = d-2,d>2; 

/r^xrI^I^I^(^‘-^°)I'’ ifd-2<s<d,s>Q; 


where both sides remain finite. Also, if p° is compactly supported, then p* remains compactly supported for 
all t > 0. 


In the case d—2 < s < d, s > 0, uniqueness remains an open problem: it has been obtained in dimension 1 by [5] 
(integrating the equation with respect to x and then considering viscosity solutions), but in higher dimensions no 
result is known (cf. [9, 8]). Nevertheless, as a consequence of the stability result of Lemma 2.2, we easily find that 
uniqueness of bounded weak solutions always follows from the existence of a smoother solution, so the problem 
is somehow reduced to a regularity question: 

2. This follows from the observation that for 0 < s < d — 2 we have Ags{x) = —s(d — 2 — , and hence for all p > 1 

we find (formally) along solutions 

dt I{y^Y = {p-l) j (/)PAh* < 0 . 
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Corollary 2.4 (Weak-strong uniqueness for the fractional porous medium equation). Let 0 V (c? — 2) < s < d. 
Let fi° G P(R'^) n L°°(]R‘^), and in the case s = 0 also assume log(2 -|- \x\)\fj,°{x)\dx < oo. Assume that 
equation (1.3) admits a weak solution p* with initial condition jjL, such that ii,'V^h G L°°([0,T];L°°(R'^)) for 
some T > 0. Then, p is the unique weak solution to (1.3) up to time T in the class L°°([0,T];L'^)®'^)). 

Proof. Let pf be a weak solution to (1.3) as in the statement, and let denote another weak solution to (1.3) in 
L°°([0, T]; L°°(M‘^)). By Lemma 2.2, we may then conclude 

[ [ gsix - y)d{fi* - iy*){x)d{n* - iy^){y) < 0, (2.9) 

JR‘‘ JR<i 


for all t G [0,T]. For d — 2 < s < d, s > 0, Proposition 2.3(ii) gives V 5 s*(p* — 1 /‘) G L^(IR‘^, \^pdxdf,), so that (2.9) 
becomes by integration by parts 


i 


xM 


|^nVg«*(p‘-r/‘)p<0. 


This proves Vgs * y* = S/gs*v^, and hence, applying the operator — div(|,^p-) to both sides, p* = for all 
t G [0, T], We may argue similarly in the Coulomb case s = d — 2, d > 2. □ 


As the following lemma shows, the required boundedness of is implied by a sufficient amount of Holder 
regularity for yf. 


Lemma 2.5. Let OV {d — 2) < s < d. Let y G 7^(IR.‘^) fl C'°’(IR‘^) for some a > 2 — d + s, and denote by h^ := gs* y 
the associated Riesz potential (2.5) on If s = d — 1, we further assume V/i G L^‘’(R'^) for some pq < oo. 
Then, we have 


||(V/^^V^/l^)||L=»<||p||L^ + ||p||c- 


( 2 . 10 ) 


Moreover, if s = d — 2 > 0 we have ||V^h^||LP IImIIlp for all 1 < p < oo, and if s = d — 1 we have 
II <p IIVpIIlp for all po <p < oo, p > 1. 

Proof. Without loss of generality we may assume y G as the claimed result then follows by an obvious 

approximation argument. We first prove that for any y G C)P(R‘^) the Riesz potential satisfies (2.10). We only 
argue for the second gradient ||V^/i^||l“, the other part being similar and easier. Let y G be symmetric 

around 0, with y = 1 in i?i and y = 0 outside B 2 . If d — 2<s<d—1, decomposing 


V^hf^ix) = [ gs{x - y)V'^y{y)dy 

= gs{x-y){l-xix-y))V^y{y)dy+ gs{x - y)x{x - y)Vl{y{y) - y{x))dy, 

jRd jRd 

we find by multiple integrations by parts 

V^/i'"(a;) = / V'^gsix - y){l - x{x - y))y{y)dy + / S/'^gs{x - y)x{x - y){y{y) - y{x))dy 

- y{x) / gs{x - y)y^xix - y)dy, 

Jr'^ 

and hence, for any 2 — d-|-s<CT< 1, 

|V^/i^(x)| < /" \x - y\~"~^\y{y)\dy + \\y\\c- [ \x-y\'^-‘‘-^dy+\\y\\i^<^ [ \x-y\-^dy 

Jlx — vl^l J\x — v\<.2 J \x—y\<.2 


that is (2.10). If d — 1 < s < d, rather decomposing 

i / ^9s{x - 2/) 0 (V/i( 2 /) - Vfi{2x - y))dy 

= / ^^gsix-y){l-x{x-y))fj.{y)dy- Vgs{x - y) ®Vx{x - y)y,{y)dy 

JvA dRd 

+ \ j Vgs(a; - y) 0 g[y) - Vy.{2x - y))x{x - y)dy, 

^ dR<J 

the result (2.10) again follows from a direct computation. As far as the additional L^-boundedness is concerned, 
it is a direct consequence of the L^-boundedness of Riesz transforms for 1 < p < oo, simply noting that we have 
_ v 2(-A)-V for s = d - 2 > 0, and ~ V(-A)-i/2Vp for s = d - 1. □ 

Motivated by these considerations, we would like to prove at least that the regularity of the initial condition is 
conserved along the flow, so that the boundedness of would be ensured by the above lemma for sufficiently 
smooth initial data. In the Coulomb case s = d — 2, d > 2, any weak solution p* to (1.3) can be shown to belong 
to L°“(R+;'P(]R‘^) n C"^(IR‘^)) whenever /i° G 'P(R'^) n C‘^(M‘^) with non-integer cr > 0 (cf. [19, Theorem 1], which 
is easily rewritten in any dimension). By a similar but easier argument, the corresponding result can also be 
proven to hold in the case 0<s<d — 2,d>3. In the case d — 2<s<d—l,s>0,a recent result [30] shows 
that there exists a time T > 0 (depending on initial data) such that any weak solution p* to (1.3) belongs to 
L°°([0, T]; iJ'^(R‘^)), whenever G T’(K'^) fl iJ'^(R'^) with tr > | -|- 1. In the case d — 1 < s < d, however, even 
such a short-time regularity result is unknown. This is why we needed to add some regularity assumption in the 
statement of Theorem 1.1. 


2.4 Modulated energy and elementary properties 

Let 0 < s < d, and let /i°, y.*, y.% be as in the statement of Theorem 1.1, for some T G (0, oo). Let A > 1. 

Since is assumed to have bounded energy, and since the energy is decreasing along the flow, we find 


N 


N 


-^^NixXN, ■ ■ • ,a;)v,Ar) = sui^^ -^^dsixlN - x\^) < ^ ^3s(<Ar - x^^) < oo. 


te[o,T] 


te[o,T] 


i¥^3 


If 0 < s < d, since gg is nonnegative, this proves 


N 

77Ar := min inf jx-^ - x) ^| > 0. 

If s = 0, po changes sign and some more work is then needed: noting that by symmetry 

1 AT ^ N ^ N 


‘^j,N 


' N ■ 

a direct computation yields 

N 




= 0 , 




dt 


N N 9 ^ 




2=1 

N 


N 2 


i¥=3 


i=l 


l<Ar "4^1^ 


N 




^i,N ^j,N _ 2(A — 1) 




\^t _ rpt 12 


N 


( 2 . 11 ) 
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and hence 


1 " 

- xIj^) + c^ 1 \xIn - xInI"^) 

tG[0,T] 

1 ^ r~^ ^ 0/AT T\ 

- 7 ^ X! 9o{xIn - x^n) + X! 


i^j 




N 


< oo. 


As go{u) + cJqVT > 0 for all u, this proves that (2.11) also holds in the case s = 0. 

Next, we recall the truncation procedure introduced in [21], which serves to make energies finite without 
removing the diagonal. For fixed TV > 1, let 77 > 0 be small enough such that 277 < 1 A tjn- Then define 

N 






N 


where denotes the uniform unit Dirac mass on the sphere dB[z,r]). Denote for simplicity 


h gs * g, , hpf gs * g,pf, 


f^N,n 9s * Mtv,) 


and use the same notation for their extensions to x R as in Section 2.2. Also define := gs{n) A gs- Noting 
that by symmetry 

-divd^pVffg.,,) = 4’’^4<1x{o}) 

where <5o**^4‘ix{o} denotes the unit Dirac mass on dB^j x {0}, we find 


N 


^N,rj{x,^) — ,j. ^ ^ 9s,r){x 


N 

I—I 

Let us now introduce our notation for the small balls around the particle locations, which we will be crucially 
using in the proof: for all i? > 0, let B%{R) denote a union of disjoint balls 


Ml,(R) 

Blf{R) := [J 

m—1 


( 2 . 12 ) 


with total radius R = n such that x* pf € B%{R) for all 1 < T < A^. These balls will be carefully 

chosen in Section 2.7 below. 

As already announced, for all > 1, we will consider the following modulated energy 


IJd 9six-y)d{fi% - y*){x)d{fi% - y*){y), (2.13) 

where D := {{x,x) : x G R'^} denotes the diagonal. This quantity can be thought of as a natural renormalization 
of the iJ“*^'^“®)/^-distance in the presence of Dirac masses. Its main property is as follows: 

Lemma 2.6 (Modulated energy). For allt >0, if the sequence is tight, then the following two conditions 

are equivalent: 

(i) limsup^|^£’Ar(t) < 0; 

(ii) gf-pp ^ g* and gs(x - y)dg%{x)dg%{y) gs{x - y)dg*{x)dg\y). 

Proof. Property (ii) clearly implies (i) (and even £N{t) —>■ 0), so it suffices to check the converse. Assume that 
limsup^ fAr(t) < 0. By tightness, up to extraction of a subsequence, the Prokhorov theorem gives g*pf ^ v* for 
some v* € P(R‘^). For any AT > 0, we may write 



y)dglf{x)dg%{y) > K A gs{x - y)dg%{x)dglf{y) 

= -^+[ [ K Ags{x-y)dg%{x)dg%{y), 

tv Js.d Js_d 
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and hence, successively passing to the limits iV f oo and if f oo, we find 


9 six-y)d^J,%{x)d^l%{y)> / gs{x - y)dy\x)dv\y). (2.14) 

^too JJjjc J^d J^d 

Combining this with convergence /i)^ and with assumption limsup^ £N{t) < 0, we obtain 

0 >limsup// ge{x-y)dy,%{x)dfi%{y) -2 / gs{x-y)diy\x)dy*{y)+ / gs{x - y)dy\x)dg.*(y) 

Wtoo jRd JwLd J-Rd iRd 

>/ / gs{x-y)dty\x)di'*{y) -2 / gs{x-y)di'*{x)dy\y)+ / gs{x - y)dg*[x)dg*(jj) 

= [ [ 9s{x-y)d{v* - g^){x)d{v* - g^){y). 

JR'i jRd 

The result then follows, noting that y* has bounded energy by Proposition 2.3, that has bounded energy 
by (2.14), and noting that for any two Radon measures /i, u with hnite energy we have 

/ / 9 s{x-y)d{v - y){x)d{v - g){y)>Q, 

JRd JR'i 


with equality if only \i y = v (see e.g. [18, Theorem 9.8] for 0 < s < d, and [22, Lemma 1.8] for s = 0). 


□ 


In the case of bounded weak solutions to (1.3) as given by Proposition 2.3, the following identity follows 

from an integration by parts and was crucially used in the proofs of Lemmas 2.1 and 2.2 (cf. (2.4) and (2.8)): 


/ f- 


gsix - y)d{y\ - y\){x)d{g\ - y\){y) 


4 ,|V(h‘ -h*)| 2 , ifs = d- 2 , d> 2 ; 

/r-^xR l^ri^(M - if d - 2 < s < d, s > 0 . 


(2.15) 


Now we would need a corresponding identity in the context of the modulated energy £]si(t). Since does not 
belong to L^(K‘^) or L^(K'^ x K, |^|'’'d 2 ;d^), a regularization is then needed. Besides the modulated energy Spf, we 
thus define the following 77 -approximation, based on the truncation introduced above: 


g U-. .= J if s = d - 2 , d > 2 ; 

l/R^xRl?r|V(^kr 7 -^‘)P^ ifd- 2 <S<d, S> 0 . 

An integration by parts then yields the following proxy for identity (2.15), showing that the difference between 
the modulated energy £N(t) and its approximation £is[,r](t) just comes from the diagonal terms (which are indeed 
excluded in £Nit) but not in £N^^{t)). We refer to ]21, Section 2.1] for a detailed proof. 

Lemma 2.7 (Approximate modulated energy). Let 0 V (d — 2) < s < d. For all t > 0, N > 1 and g > 0, 

£N,r,(t) = £N(t) -f + O^^^(l), 

where for any fixed N we have o^^^(l) —> 0 as 77 j, 0. 


2.5 Gronwall argument on the modulated energy 

By Lemma 2.6, in order to prove convergence -^^011 ^g convergence of energies, up to tightness 

issues, it suffices to check that limsupjy fAr(t) < 0. This is achieved by a Gronwall argument. From now on we 
focus on the Riesz case d — 2<s<d, s>0. The Coulomb case s = d— 2, d>2 can be treated in exactly the 
same way, but is actually easier since it does not require to use the extension representation of Section 2.2. 
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Proposition 2.8. Let d — 2 < s < d, s>0. Let be a probability measure such that equation (1.3) admits a 
solution jj} that belongs to L°°([0, T]; P(R'^) nC'^(R'^)) for some T > 0 and some a > 2 — d+ s. In the logarithmic 
case s = 0, d = 1, also assume that /jjjlog(2 + \x\)\gL°{x)\dx < oo, and V/i G L°°([0,T];L^(]R)) for some p < oo. 
Let pifi ^ pL°, assume 


limsup // gs{x - y)dfj.%{x)dp,%{y) < oo, 
N'\'oo JJ x^y 


(2.16) 


and let /i)^ be defined by (1.1)-(1.2). Assume that, for all t G [OjT], the collection can be chosen with 

0 in such a way that 


lim inf lim inf 

N'^OC 774'0 


t ,2 9s{v) 




icriv/iVJ- 


N 


> 0 , 


and, denoting gf ft) '■= c^].t ® for s > 0 and g^ ft) := (,(—logt) V 0 otherwise. 


1 


N 


}}^JpJ29t{d{xlN,dB*f^{R%)))=0. 


2=1 


(2.17) 


(2.18) 


Then, for all t G [0,T], we have SNft) < Ct(^Ar(0) + OAr(l)). 

Remark 2.9. In the ideal case when all particles remain well-separated, that is with a minimal distance rjN of 
order then, taking B%{R%) to be the union of balls of radius R%/N centered at the points xl jy’s, with 

R%/N <C condition (2.18) simply becomes gs[R\q/N)/N 1. On the other hand, neglecting interactions 

between particles, hence focusing on the (divergent) self-interactions, we formally find 





N 

^ ' -- — I 'I \x — x. 


1 

iv 

1 

'n 


i=l 

[ i^riv5.(x,oi^ 

Jri<\x\<R%/N 

{9s{v) ~ 9s{RN/dI)) + •. •, 




so that condition (2.17) would amount to requiring gs{R%/N)/N 1, which is thus just the same as condi¬ 
tion (2.18). In other words, for s > 0, both conditions would then take the form ^ which is 

compatible with R)^ —)> 0 only if s < 1. In Section 2.7, we prove that a consistent choice of the small balls 
B*N{Rf^) is indeed possible whenever 0 < s < 1. 

To go beyond the restriction s < 1 via this approach, we would need to modify Proposition 2.8, in particular 
by refining the (blind) approximation argument used in Step 2 of the proof below, in order to relax the smallness 
condition for the total radius Rf^ —>■ 0. To do that, precise microscopic information on the particle dynamics 
would become needed. Getting a handle on such information seems however to be a difficult task and is not 
pursued here. 

Proof. By the regularity assumption for /i*. Lemma 2.5 ensures that we have ||(Vx/i*, V^/i‘)||l=“ St 1; and also, 
in the case s = 0, d = 1, ||V^/i*||lp St 1 for some p < oo. We split the proof into four steps. 

Step 1: Time-derivative of S^ft) and modulated stress-energy tensor. In this step, we prove equality 


dt£N{t) = - [ m^VlhSx) : TS{x,0dxdf 

Jr‘‘xR 

-2 f p.v. f VgSx-y)d{p.% - p!^)fy) 


dpStix), 


(2.19) 
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:= lim 

x} ’■'I'O jR‘i\B{x,r) 


where we use the usual principal value symbol 

P.V./ 

and where the modulated stress-energy tensor is defined as follows: for all 1 < fc, 1 < d, 

■= ^ JJ dkgs(.x-y,^)digs{x-z,i)d{^iN - ^J){y)d{^iN - ^J){z) ( 2 . 20 ) 

- 4; // ygs{x - 2/, 0 • '^gs{x - z, ^)d{g,N - g){y)d{g.N - g-){z). 


£)■= 

Moreover, as checked at the end of this step, the integrals in (2.19) are summable: more precisely, we prove that 
|r)(f| belongs to L^(K'^ x R, \^\zdxd^) if s > 0, and that |V^/i*(x)||r^(a:,^)| belongs to L^(R'^ x R, \^\zdxd^) if 
s = 0, d = 1. Although the second term in the right-hand side of (2.19) is nonpositive, we do not bound it by 0 
yet, contrarily to what is done in the proof of Lemmas 2.1 and 2.2, since it will be useful in Step 2 below to absorb 
some error terms. 

Using the equations satisfied by /i‘ and by the trajectories xl and noting that the gradient Vh* is given by 

Vh\x) = p. V. f ygs{x - y)dg.\y), 

Jv.'^\{x} 

where the principal value may only be omitted for s < d — 1, we find the following expression for the time-derivative 
of the modulated energy £’Ar(t) defined in (2.13): 

1 2 ^ 

dtSN{t) =dt / gs{x- y)dg* {x)dg* {y) + dt^^gs{x\N - x] - 4^ V / gsix^N - y)dg^{y) 
jRd dRrf ^ A ^ djjd 

l^J 1^1 

N 

\[ yh\x)-p.Y. [ ygs{x-y)dy\y)dg*[x)-^^ ^^Vgs{xljk,-x]j^) 


= -2 


2 ^ 2 ^ 
^ -Vgsixl^N - x]^n) + 




^\.n} 


Let us rearrange the terms as follows: 

dt£Ni.t) = -2[ p.v. f Vgs{x-y)d{fi% - g.*){y) 


Vh\x) ■ Vgs{x - xl^)dg,\x). 

I 

dy%{x) 


-2 

-f 2 

-2 

-f 2 


[ Vh\x)-p.v.f Vgs{x - y)dg,\y)dy*{x) 

Js.'^ iR'^\{a:} 

Vh*{x)-p.v. / Vgsix-y)dg.\y)dy%{x) 

Js.'‘\{x} 

Vh*{x)- Vgs{x-y)dy%{y)dg.%{x) 

Jr^\{x\ 


R'i 


p.v. 


/R‘1\{x} 

V/i*(x) • Vgsix - y)fi*{x)n%{y), 


/R<i\{y} 

and note that the last four terms in the right-hand side may be combined to yield the following simpler expression: 

2 

dg,%{x) ( 2 . 21 ) 


dt£N{t) =-2 f p.v. [ ygs{x - y)d{yL% - y^){y) 

dRd Js.<^\{x} 

- “ Vd*( 2 /)) • Vgs{x - y)d{y% - y*){y)d{^i% - g^){x) . 


— : lN{t) 
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In the distributional sense on using canonical regularizations, we may alternatively write 


lN{t) = {S%;Vh^) = dkh*), ( 2 . 22 ) 

k^l 

where and, for all 1 < /c < d, 

[ dk9s{x-y)d{^i%- 

JR'‘\{x} 

Since — div(|^pV 5 s(x —xq, ^)) = (5a,o(3^)<5B‘^x{o}(2^) 0 for xq € M'^, we have in the distributional sense on xR 
x{o}(a^> 0 = -2 P- V. div(|^|T'V5s(a; - z, 0)9k9s{x - y, - y*){z)d{y% - y*){y) 

II “ z,0)dk9s{x - y,0 +div{\^pVgs{x - y,^))dk9s{x - z,^)) 

X d{y% - y*){z)d{y% - 9*){y). 

Now note the following algebraic identity in the distributional sense on R'^ x R: for all 1 < /c < d, 
div(|^pVgs(a; - y,Cj)9k9s{x -z,^) + div(|^pV5s(a; - z,^))dk9s{x - y,0 

1 d-\-l 

where we have set 

tZ+l 

Gi^{x,^;y,z) := 2di9s{x - y,^)dk9s{x - z,^)- Sik ^ djn9s{x - y,0dm9s{x - z,^). (2.23) 

m—1 

This proves the (Delort-type) identity 

d-\-l 

Sf{x)6^<^^mix,0 = - (2.24) 

1=1 

for all 1 < fc < d, and the conclusion (2.19) then follows from (2.21), (2.22) and an integration by parts. 

We now turn to the claimed integrability of the modulated stress-energy tensor T^. We first consider the case 
d — 2<s<d, s>0. For that purpose, we begin with the bound 

[ ^ [ \^P [[ \{x-y,0\~"~^\ix-z,^)\-^-^d{9%+fi^){y)d{n%+f/){z)dxd^ 

dR>ixR dR-^xR JJ D'= 

= [f ([ \^P\{x-y,0\~'‘~^\{x-z,^)\~'‘~^dxd^'\d{9%+y^){y)d{n%+f/){z). (2.25) 

JJD'= VdR<ixR / 

Let us compute the integral over R'^ x R. Denoting for simplicity Cyz := {y + z)/2 and g := s -I- 1, we decompose, 
for all y z, 

[ (|a; - dP + - z\^ + \)-^Gdx = -p -p -p 

dR<i 
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where 


/I ■= 

yz 


r := 


r ■= 

yz 


:= 
yz 


'k-y|<||y-2| 

:= [ {\x - yp + l)-^/"(|x - zp + 

i|a;-2|<i|y-z| 

X - yp + l)-«/2(|x - zp + 

(|x - yp + l)-«/2(|x - zp + l)-«/2dx, 


-I- 

< 


-y|,|x-z|>^ |y-z| 
\^-cyz\<\y-z\ 


’\x-Cy:,\>\y-z\ 

Using that |x — y| < ^|y — z| implies |x — z| > ||y — z|, we may estimate 


4 < (12/ - ^1V4 + 1)-^/" [ {\x - yp + < (|y - z|/2 + 

y |x-y|<i|y-z| 

and similarly for ly^. Moreover, 


(|y - z|V4 + l)-« [ dx< (|y - z\^/A + l)-^|y - z\‘^ < (|y - z|/2 + 


lyz < ( 12 / - 

'\x-Cy^\<\y-z\ 

and also, since d — 2q < 0 follows from the choice s > d — 2, s > 0, 


r < 

yz 


< 


\x-Cy^\>\y-z\ 


\x-Cy^\>\y-z\ 


(|x — y| + 1) ‘^{\x — z| + 1) ’^dx 
I 

(|a: - Cy^l - |y - z|/2 + l)~'^’^dx < (|y - z|/2 + 


This proves, for all y ^2 2 ;, 

[ (|x - yp + l)-«/"(|x - zp + < (|y - z|/2 + 

yR'i 

and hence by scaling 

[ \{x - y,Or^|(x - z,Ordx < (|y - z|/2 + 

so that we obtain, as by definition 7 = y — ds (—1,1), 

[ \^r\ix-y,0\-‘‘\{x-z,0\-‘^dxd^< [ l^l^-^^dy-^l + ieD'^-^^d^, 

yR'^xR yR 


Splitting the integrals over ^ into the part where |^| < |y — z:] and that where |^| > |y — z|, and noting that y > 1 
follows from s > 0, we find 

[ \^r\{x-y,or\i^-z,ordxd^ 

yR'^xR 

|y - [ |er-"de + [ < |y - z4« = |y - z|-^ 

'^\£\<\v—z\ ^\£\>\y—z\ 


< b/ — r., 

'lCI<|y-z| 

Combining this with (2.25) finally yields 


[ ICri4f| ^ [[ \y ~ ^\ ’^d{y,% + fi*){y)d{y% + f/){z), 
yR'^xR jj D‘^ 
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and hence, by assumption (2.16), since both the particle and the mean-field energies are decreasing along the flow 
(see Proposition 2.3 for the mean-field energy). 


[ [ [ gsiy - z)df/{y)dfi\z) + [[ gs{y - z)dy%{y)dy%{z) + 2 f h*dy,% 

jR'ixK JR'i Js.‘‘ JJD'= J 

<[ [ gsiy - z)dfi°{y)dy°{z)+ [[ g^iy - z)dy%{y)dy%{z) + 2\\h*\\i,^ <t 1. 

jRd' JRd' JJ 

We now briefly consider the case s = 0, d = 1 (hence 7 = 0, g = 1). Let 1 < p < 00 be such that 1- 

Arguing as above, we obtain 

[ iV^h^xMx - 2 /,e)|-'|(x - z,0\-^dx < ||V 2 h*||L»(| 2 / - Z| + |?|)-\ 

Jr 


and similarly, by the Holder inequality, for 1 /p -|- 1 jp' = 1, p' > 1, 


[ \v^h\xm^-y,om^-z,o\-^dx<\\s/^h*\\^J [ \{x-y,or'\{^-z,or'dx] 

dR VdR / 


1/p' 


<||v2h*||L.(|y-^l + iei)^-^ 

Splitting the integral over ^ into the part where |.f| < |p — 2 ;| V 1 and that where |^| > |p — 2 ;| V 1, we may then 
estimate 


[ \W‘^h^{x)\\{x -y,0\ ^ 1 ( 2 ;- 2 :,61 ^dxd^ 

Jrxr 


X 

2ut 


< \\V^h%,^ [ (|y_^| + |^|)-lde+||V 26 ||LP 

kl<|y-2|vi 

<t 1 - 0 A log(|p - z|) = 1 -h 0 V go{y - z), 


i 


|C|>|y-z|Vl 




so that the conclusion now easily follows just as in the case s > 0 . 


Step 2: Approximation argument. For all t > 0 and all R S (0,1), applying [25, Proposition 9.6], there exists a 
smooth approximation u* of the function V/i* G (7°’^ (R"^; R'^) such that v* is constant on each ball of the collection 
satisfies, for all a G [ 0 , 1 ], 

||6 - V 6 ||c- < C'i?i"“||V 2 /i‘||L- < CtR^-°‘, (2.26) 

and also satisfies ||Vu*||lp 1 for some p < cxd in the case s = 0, d = 1. In this step, we prove the following 

estimate: 


dtSNit) < - [ + Ct or{1), (2.27) 

dR^xR 

where 0 ^( 1 ) denotes a quantity that goes to 0 as i? 0 . 

Using relation (2.24) as well as the integrability properties of T^, we may decompose the first term in the 
right-hand side of (2.19) as follows: 


[ (C(6V.6(x)) : r^(x, Odxd^ = {s%-, V6) = {s%-, u‘) + {S%; V6 - u‘) 

dR<lxR 

= {^n^R'‘x{o}tX*) + {S%; Vhd — V*) 

= f :r^ + (56;V6-r;‘). (2.28) 

dR'^xR 
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It remains to estimate the last term in the right-hand side of (2.28). Denoting for simplicity uf := V/i* — t;*, we 
may decompose 


(2.29) 


{Si,- - v^) = - w\y)) ■ Vgs{x - y)d{gl - y*){y)d{y*,, - y^){x) 

= 2/ [ w*{x)-Vgsix-y)df/{y)dy\x) 

Js.'i' Jr‘‘ 

-4 [ [ {w\x)-w\y))-Wgs{x-y)dg%{y)dg\x) 

Js.'i JRd 

+ 2^ w\x) ■ Vgs{x - y)dij,l,{y)dij.l,{x). 

For the first term in the right-hand side of (2.29), we simply have by (2.26) 

[ [ w*{x)-Vgsix-y)df/{y)dg\x) = [ w*-Vh^df/ < ||w‘||l“ || V/i‘||l- < Cti?. 

jRd' JR'i JRrf 

As far as the second term is concerned, choosing (t>s-|- 1 — d, 0<cr< 1, and recalling that remains bounded 
by assumption, we find by (2.26) 


/ / {x) - (y)) ■ Vgs(x - y)dgl{y)dg^(x) 

JR‘‘ JRd- 

I \x-yr-^+'^dy\x) 


^ \\w lie- sup 

1/gR'i 

< ||rc*||c- sup 


J la 


|x-y|<l 


\x - y\ 


— S —l+(7 


dx 


[ dg\x)] 

d|x-y|>l / 


< llw^lle-)! + ||d*||L°°) < CtR^ ■ 

Combining these two bounds with (2.29), using (2.26) once again, we obtain, for i? J, 0, 


I (5: 


-v*)\<tOR{l)+R [ [ Vgs{x-y)dgliy) 

Jr'* Jr‘^\{x} 


dgl{x) 




or{1)+r[ p.v. / Vgs{x-y)d{g% - g.*){y) 
Jr<‘ Jr<‘\{x} 


dg%{x) 


R 


R<i 


p.v. / Vgs{x-y)dg\y) 


dgl{x), 


and thus, noting that 


/ P-v. / ygs{x-y)dg\y) 
dRd dR<i\{x} 


dMx) < ||Vd‘||L~ <Ct, 


we find 


i(^: 


-v*')\<tOR{l)+R [ p.v. [ Vgs{x-y)d{g% - g*){y) 


/R‘i\{x} 

Hence, for all e G (0,1), using inequality i?|a| < ea^ -I- (4e)“^i?^, we obtain 

\{Sl-,Vh^ -v*)\S<^~^OR{l) + e f p.v. [ y gs{x - y)d{g% - g^){y) 

dRrf JR<i\{x} 


dgl{x). 


dg%{x), 


and the result (2.27) then follows from (2.28) and (2.19), choosing e > 0 small enough (depending on t). 
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Step 3: Modification with rj-approximations. In the definition (2.20) of T^, the diagonal terms were excluded. In 
order to apply inequality 2|a6| < to as in the proof of Lemmas 2.1 and 2.2, we would need to add 

these diagonal terms explicitly. Moreover, 77 -approximations then become needed in order to avoid the divergence 
of the corresponding diagonal terms that will appear after application of the above-mentioned inequality. More 
precisely, we prove in this step 


dt£N{t) <t 


1 




(2.30) 




|^|^|V7;‘(x)||V5s(a:: - + Oij(l) -f 


J {R'^\B]^{R))xR 

By the choice of v* to be constant on each ball of the collection B\^{R) and the bound on Vu*, (2.27) becomes 

dt£N{t)<t[ |enV7;*(x)||r^|+Ofl(l). (2.31) 

J (R<‘\Bi,(R))xR 


Denote for simplicity 


H%ix, 0 := - h*){x, 0, Jx, f) := {h% - h*){x, 0, 


and define, for all 1 < fc < d -|- 1 and all 1 < Z < d, = 0 and 


■■= 2dkHN,^{x,f)diHN,^ix,0 - 


For all X with d{x, > 77 , we note that 

N 

N 


N 


= ^'^Vgs,^{x-xlN)-Vh\x) = ^J2'^9s{x-xl^)-Wh\x) = VH%{x,f). (2.32) 


Also noting that definition (2.23) may be rewritten as 


1 

iv2 


xl^,xl^) = 2 [[ dk9s{x - y] f)digs{x - z-0d{tJ*N - h^){y)d(.p% “ 
i=i jjd 


- Ski II ^gsix - y; 5) ■ ^9six - z; - n*){y)d{p% - y*)(z). 


D 


definition (2.20) yields 

1 ^ 


= 2 


/ / dkgsix-y;09igs{x - z;C)d{y% - y,*){y)d{y.% - g*){z) 

JR'^’ dR<i 

-Ski / ygs{x-y;£,)-"^gsix - z;^)d{f4r - y‘){y)d{y% - y*){z) 

Jr'^ Jr'^ 

= 2dkH%{x,f)diH%{x,f) - Ski\'^H%{x,f)\'^. 

Combining this with (2.32) yields, for all 1 < Zc < d + 1, all 1 < Z < d, and all x with d{x, {x* jvlfci) > V-i 

1 ^ 


2=1 
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From (2.31), we then deduce, for all 77 > 0 small enough such that (Jili v) C B%{R), 


N 

dtSN{t)S f mTh.r,\ + ^^[ \^r\Vv\x)\\Gs{x,^;xy,xl^)\+OR{l). 

JiR'‘\B%{R))xR d(Rd\S^(iJ))xR 

The result (2.30) then follows, using inequality 2|a6| < + 6^ in the form of 

I^AT.r,! ^ - h*)p, \Gsix,^;xl^,xlj^)\ < |Vg«(a; - OP- 

Step 4- Conclusion. In this step, we show that 

dt£N{t) S £Nit) +(^- [ 

\ dB‘ (fl)xR 


•nW' 

N 


(2.33) 


1 

+ ■^J29tid{<,N,dB%{R))) + or{1) + of’«)(!). 


The statement of Proposition 2.8 immediately follows from this inequality, with the suitable choice of i? = i?)^, 
together with a simple Gronwall argument. 

By Lemma 2.7, inequality (2.30) may be rewritten as follows: 

9s{v) _ 

/B^(it)xR 

+—'yf 

^ i(R‘i\B*,(ii))xR 

or equivalently, expanding the square. 


dt£N{t) s £N{t) + ^ / l?riv(hf ^ - h‘)p 


|Cr|Vu‘(a;)||Vg,(a; - f ^,C)pda:d5 + or(1) + of’-^^(l). 


dt£N{t) <t £N{t) + d-M. 


[ i^rivhf^p-/ mvhr+2[ 

JB%{R)x'R. Jb%{R)x'R. Jb%{R)x'R. 


N 


+ 7V2 ^ 


J{R‘^\B%iR))xR 


|^P|Vz;‘(x)||Vg^(x - xli^,^)\'^dxd^ + Oii(l) + of’^^(1). (2.34) 


The last term in the first line of (2.34) is easily estimated as follows, using the notation (2.12) for the union 

/ -r~^\ r 1 


B]^{R) of small balls. 


[ i^rvhf,, • vh‘ 

7 s* (R)xR 


< l|Vh*| 


N 

l?ri(®-fAr>OI ^dxd^ 

^ :^7b*^(r)xR 


N 


<. 


-V 

N ^ 


i=i Jb*n(R) 


lx — ‘^dx 


Ml 


Ml 


\xYdx < y 

while the term in the second line of (2.34) is, in the case s > 0, 


Xjn.N dif 


(2.35) 


m—1 


1 f 

/ l?r|Vu‘(x)||V5„(x-xy,C)pdxdC 

^ i=l 7 (R‘*\B^(fl))xR 




< 


1 ^ 


iV2 


i=l 7(Rd\B^(fl))xR 
N 




X 


7r<*\B*,(R) 


1 ^ 

-xyi "dx < — y d(xy,ai3)v(i?))’ 


(2.36) 


2 = 1 
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In the case s = 0, d = 1 (so 7 = 0), denoting p\ ^ := d{x\ ^, dB*'j^{R)), and applying the Holder inequality with 
1/p + 1/p' = 1, where p < 00 is chosen in such a way that || Vu^IIl? 1, the term in the second line of (2.34) is 


N 

Wa 


< 


< 




1 

JP 

1 

Jp 

1 

Jp 


J(R\B%iR))xR 
N „ 


|Vu‘(x)||Vpo(a:: - xl ff,^)\'^dxd^ 


E 


N 


\Vv\. 


x}\\x — X. 


i,N\ 


^dx 


E/ \'^v\x)\\x-xIn\ + \Vv\x)\\x-xIn\ ^dx 

N 

u 


i=l ■'Pt,N<\x:-xltf\P 

and hence, by the choice p' > 1 , 

jpi 


\x-x\^\ dx + 




= 1 V „|>1 


\x - x\ j^\~^ dx] , 


\Vv*{x)\\\'goix-xljs!,0\'^dxd^<t ^^i-0Alogpl^)+N (2.37) 


N 


7V2 


{R\Bj^{R))xS. -■ 

The result (2.33) then follows from inequality (2.34) together with (2.35) and with (2.36) or (2.37). □ 


2.6 Bypass of tightness issues 

Assuming that fAr(O) < OAr(l), Proposition 2.8 yields £N{t) < C'tOAr(l) for all t € [0,T]. If we know that the 
sequence {p%)n is tight, Lemma 2.6 then allows us to conclude with the desired convergence ^ /i*, while 
tightness can easily be checked under the additional assumption that the initial measures p^’s are well localized 
in the sense that limsup^ J |a;pdp^ < 00. However, in the spirit of [28, Section 4.3.5], the following refinement 
of Lemma 2.6 shows that much more information may be directly extracted from the fact that £N{t) < C'tOAr(l), 
so that in particular tightness is obtained a posteriori without any additional assumption. 

Corollary 2.10. Let the assumptions of Proposition 2.8 prevail. Also assume that £n(0) < Oiv(l)- Then, for all 
t G [0,T], we have —>■ Vh* in L]/^(R'^; L^(IR, j^l^d^)) for all 1 < p < 2d/{s + d), and hence ^ p*. In 

particular, (p)y^)jv is tight and Lemma 2.6 then implies the convergence of the energy. 

Proof. By assumption. Proposition 2.8 yields £N(t) on{I). We split the proof into three steps. 

Step 1: Strong convergence outside small balls. In this step, we prove 

[[ l?r|V(h‘v-h‘)pSoiv(l), (2.38) 

dd(R<i\B^)xR 

and hence, for any 1 < p < 2, the Holder inequality implies for all i? > 0 

f (f ieriv(h)v - ([[ icnv(h)v - <«.* o^(i). (2.39) 

Jbr\b*^ \Jr J \JJ{rPb*^)xr / 

Applying Lemma 2.7 and expanding the square, we may decompose as follows the L^(]R‘^ x K, |^|^da;d5)-norm 
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of V(/i^ 17 ~ outside the small balls 




SnA^) - // 

ieriv(h* -h‘)|2 





£»(f) + A 

1 

X 

ienv(h)Y,,-/i‘)p 


«»(*) + A 

-// 

JJB%xR 

ierivh)Y.,p- [[ 
JJK 

i^rivh‘ 

, xR 

.Nf 

JjBirXR 

lervh^v,. 

■ V/r‘ + oW(l). 



Applying Proposition 2.8 in the form of Snit) OAr(l), and using assumption (2.17), this turns into 
[[ l?r|V(/i)v., - h^)\" < 2 [[ |crvh*v,, • V/i* + CtONil) + oW(l). 

dd(R<i\B^)xR JJBlfXS. 

Now arguing just as in (2.35), we find 

[[ • vh* 

-'de^xR 


^7V OAr(l), 


and hence 


// 

JJ a 


|^|7|V(/,* -/7*)|2 <, 07^(1) +oW(l). 


(R-^VZS^ixR 

Passing to the limit 77 j, 0 in this inequality, and noting that “t V/i)y hr the distributional sense, the 

result (2.38) follows. 

Step 2: Neglecting the contribution inside small balls. The contribution inside the small balls B% is of course infi¬ 
nite, since does not belong to L^(R'^xR, \^\'^dxd^). However, we show that it is small in L^(R, |^|'’'(i'f)) 

for p small enough. More precisely, for any 1 < p < 2d/{s + d), we show that we have for all i? > 0 


JBRnB* \Js. / 


(2.40) 


Decomposing Vh*j^{x) = ^ X]i=i — xl the triangle inequality yields 

p/2xl/p ^ N 


\ p/2 \ 1/p 




A direct computation of the integral over ^ yields 

P/2\ 1/p 




leriV/r^v'^ 


< 


N 


jzif \X - 

i=i \Jbr J 


1/p 


As for each i the integral over x € is clearly bounded above by the same integral over x S B{i{x\ ^y), we obtain 


IbrhB* 


[ krivh 

Js. 


t |2 

n\ 


p/2\ 1/p 


< 


(/ 

\Jbr 


|f(7+l-2(s+l))^^ 


1/p 
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and hence, for any 1 < p < 2d/{s + d), 

Jerdb^ \Js. ) Jbr Jbr 

Now, for any 1 < p < 2d/{s + d), choosing any p < q < 2d/{s + d), the Holder inequality yields 
r f f \P/2 / f / r \9/2\p/g 

/ / i^rivh^viM / / i^nvh*viM 

JBRnBl, \Js. J KJBRnB]^ \Jr J J 

The result (2.40) follows from this and from the Holder inequality in the form 

/ r \ p/2 / /’/’ \ 

/ i^nvh^n // lerivh^n <^.* 0 ^( 1 )- 

\Jr j VdiSRXR / 

Step 3: Conclusion. 

Combining (2.39) and (2.40), for any 1 < p < 2d/{s + d), we may conclude, for all i? > 0, 

[ (f ieriv(dv-d‘)p)"^'<^,,o^(i). 

Jbr \Jr / 

This proves —>■ Vd* in L^(R, |^|'’'d^)) for any 1 < p < 2d/(s + d). Applying the operator — div(|^p-) 

to both sides, we deduce p* in the distributional sense on x R, so the result is proven. □ 


L 


BRnB* 


2.7 Ball construction 


In this section, we make the heuristics of Remark 2.9 rigorous, showing that for 0 < s < 1 the collection 
can indeed be chosen with 0 in such a way that both conditions (2.17) and (2.18) are satisfied. 

Let us first describe the construction that we will use for the collection H^(i?), i? > 0. This is precisely the 
same construction as the one used e.g. in [25, Chapter 4], which was introduced by [24, 16[ for the analysis of 
the Ginzburg-Landau vortices. We first consider N disjoint small balls centered at the points x* ^’s with equal 
radii (smaller than pjv/2), and we grow their radii by the same multiplicative factor. At some point during this 
growth process, two (or more) balls may become tangent to one another. We then merge them into a bigger 
ball: if tangent balls are of the form B{ai, r^), we merge them into aiVi/ r^, ri). If the resulting ball 

intersects other balls, we proceed to another merging, and so on, until all the balls are again disjoint. Then again 
we grow all the resulting radii by a multiplicative factor, etc., and we stop when the total radius R is the one 
desired. 

As we will see, condition (2.18) is easily checked directly from the construction above, so we may focus on the 
validity of condition (2.17). Hence we need to study integrals of the form Jg* (i{)xR i? > 0. For 

that purpose, the basic tool is then the following crucial lower bound, which is a refinement of [21, Lemma 2.2[. 
In the sequel, for x G R'^ and t > 0, we denote by B'{x, t) the ball of radius t centered at (x, 0) in R'^ x R, and 
we set B[ := R'(0,t). 

Lemma 2.11 (Embryo of a lower bound). Let R > r, let be a collection of points inside the ball Br, and 

let {zk+i)l^i be a collection of points outside the ball Br. Then 


’B'r\B', 


ler 


k+l 

E 

2=1 


Vps(x - Zi,£,) 


dxdf > k^{gs{r) - gs{R)). 


(2.41) 


The same remains true if point charges are smeared out on small spheres around them, i.e. if gs is replaced by 
gs.r) with 7j < d({zi}^^l, Br \ Br). In particular, for any Zi, any R > rj > 0, and any collection (zi)j ^^2 of points 
outside the ball B{zi, R + rj). 


Ib'(zi,R) 


i^r 


1+/ 

'^'^gs.vix - Zi,f) 
2=1 


2 

dxdf > gsijg) 


9s{R). 


(2.42) 
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Proof. Step 1: Explicit value of ca.s- We claim that the normalization constant Cd,s for the Riesz kernel ps is given 
by the following formula, in terms of the beta function B(a,5) = r(a)r(6)/r(a + h) and of the measure uid-i of 
the unit sphere of dimension d — 1, 


^d,s — ^^d—1 R 


s -f 2 — d d 
2 ’ 2 


(2.43) 


Integrating the equality — div(|^|'''Vgs) = <^o on the infinite cylinder Cq := Bi x M. in x R, we find by 
integration by parts 


-f = 


[ divd^pV^s) = / \^rri-Vgs= [ f \^\'^drgs{u,Odcr{u)d^ 

J Oq QCq — oo J dB\ 

poo 

= 2uJd-i Cdrgsil,OdC 
Jo 

Since by definition gs{x,^) = cf\{\x\'^ + computing the radial derivative yields 

pCO pOO 

cd,s = 2sodd-i / r(i+ er^'^-^d^ = 

JO 

The result (2.43) then easily follows using the formula B(a, h) = /p°° f““^(l + t)~°'~^dt for all a, 6 > 0. 

Step 2: Conclusion. 

For all t > 0, denote by the ball of radius t in R"^ x R, and also fik,i ■= ^zi- We may then estimate 
by the Cauchy-Schwarz inequality 


f \^PN9s*dk,i\^ = r dt [ 

Jb'^\B’^ Jr JdBi 

>rdtf[ 1^1") Y /1^1 

Jr \JdB' J \JdB', 


'n ■ Vps * gk,i , 


where an integration by parts yields, for all r < f < i?, 

[ \^Pn-Vgs*g.k,i= [ div{\^pVgs * fik.i) =-dk,i{Bt) =-k, 
JdBi J Bi 


and hence 


f \^r\vgs*dk,if>k^rdt(f . 

Jb'\B' Jr \JdBi J 


f\^\ 

In order to compute this last integral, we use spherical coordinates: 


[ If r = [ {sm9)‘^ cos6\'^d6 = C~^^uid-i [ (1 — 

JdB[ Jo J-l 


= C^^ujd-i 


fo-ur- 

Jo 


■2)/2^(7-1)/2^^ = B 


s -f 2 — d d 


where the last equality follows from the formula B(a, b) = ^(1 — t)^ ^dt for all a, 6 > 0. By Step 1, this last 

expression is nothing but P'^^Cd,s/s, so that we may conclude 




Cd,s Jr 


R 


r^-^dt = k 


,r-^ -R- 


Cd.i 


= k^{9s{r) -gs{R)). 


□ 
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With this result at hand, as in [25, Chapter 4], we may now deduce the following lower bound estimate for 
the energy on the balls of the collection For logarithmic interactions (thus in particular for the Ginzburg- 

Landau vortices, as treated in [25, Chapter 4[), a particularly simple additive structure shows up, simplifying 
computations a lot; here we show that the same result still holds for all s < 1. 

Proposition 2.12 (Lower bound). Let i? > 0. If s <1, then, for all < rj < rjN f\ (R/N), 

[ > ^{gs{v)-9s{R/N)). (2.44) 

Jb%{R)xR 

Proof. We prove that, for all i? > 0, if B{y,r) is a ball belonging to the collection B].^{R) and contains n of the 
particles x\ ^’s, then 

[ > ^{9s(.v)-9Am))- (2.45) 

The desired result (2.44) indeed follows from summing the corresponding inequalities (2.45) associated with each 
ball B{y,r) of the collection B],^{R), noting that B'{y,r) C B{y,r) x M. We prove (2.45) by induction: we first 
show that it holds when B{y, r) contains only one particle cc* jy, and then that it is preserved through the growth 
process. 

First, suppose that B{y, r) is a ball of B\^{R) and contains only one particle x\ By definition we must have 
B{y, r) = B{x\ r) and cc* ^ ^ B{y, r + rf) for all j ^ i. Lemma 2.11 in the form of (2.42) then yields 

JB'iv.r) 

This proves (2.45) when B{y, r) contains only one particle x* ^y. 

Now we need to prove that (2.45) is preserved by the growth process, i.e. that it remains true through both 
expansion and merging of balls. On the one hand, suppose that, for some R > 0, B{y, r) is a ball of B]^{R) for 
which (2.45) holds, and suppose that B{y,r) inflates into B{y,ar) without merging, when passing from B].^{R) 
to B\j(aR), for some a > 1. Let n denote the number of particles in B{y,r). By definition, B{y,ar) contains 
the same number of particles, and the choice of y small enough ensures that no particle may lie in the annulus 
B{y,ar + rj) \ B{y,ar). Hence, Lemma 2.11 in the form of (2.41) yields 


[ \fP\yh%j>[ \^nvh%j^+[ ienvh)y,,i 

J B'{y,Oir) J B'(y,r) J B'(y,Oir)\B'(y,r) 


IB'{y,r) 


B’(y,ar)\B’{y,r) 


> ■^i9s{v) - 9 s{R/N)) + ^{9s{r) - gs{ar)). 


Since by definition r = nR/N, we hnd, by the choice s < 1, with gs{R/N) — gs{aR/N) > 0, 


[ ier|Vh)y,,P > ^{gsiv) - 9 s{R/N)) + ^ig^{R/N)-g,{aR/N)) 

JB'{y,ar) 

Tl Tl 

> j;^i9s{v) - gs{R/N)) + —{g,{R/N) - gs{aR/N)) 

TL 

- Jf2^9s{v) -gs{aR/N)), 

so that B{y,ar) also satisfies (2.45). 

On the other hand, suppose that B{jji,ri), i = 1,... ,k, are k disjoint balls of B%{R~) for some i? > 0, suppose 
that each of them satisfies (2.45), and suppose that these balls are merged by the growth process into a larger 
ball B{y, r), which is then disjoint of all other balls of the collection S^(i?). Denoting by Ui the number of points 
in B(yi,ri), we then find 


'B'(y,r) 




k 

2=1 ’ 


B'{yi,ri) 


l^r|Vh*y,j2> 




2=1 


Bi ]igsiv) - gsiR/N)), 


so that B{y,r) also satisfies (2.45). This completes the proof. 


□ 
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We are now in position to prove that both conditions (2.17) and (2.18) may be satisfied whenever s < 1, thus 
finishing the proof of Theorem 1.1. 

Corollary 2.13 (Checking conditions (2.17) and (2.18)). //0 < s < 1 and if is constructed as above, then 

conditions (2.17) and (2.18) are automatically satisfied for any choice 1 */0 < s < 1, and 

for any choice I if s = 0, 

Proof. On the one hand, Proposition 2.12 gives 

lim ( [ kr|V/i)v,,,P - -^9s{v)) > -^9s{R/N). 

»?4.o \Jb%{R)xR / f'' 

On the other hand, since by definition Uti at) R/^) C B%{R), we may estimate 

f ^ f 

■j^Y.9tid{xlN,d^%iR))) < j^9tiR/N). 

i=l 

Therefore, both conditions (2.17) and (2.18) are satisfied if we choose such that ;^5j~(7?)v/iV) <C 1, and the 
result follows. □ 
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